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Abstract. A new high precision measurement of the reaction pp ^ pK+A at a beam momentum of 
2.95 GeV/c with more than 200 000 analyzed events allows a detailed analysis of differential observables 
and their inter-dependencies. Correlations of the angular distributions with momenta are examined. The 
invariant mass distributions are compared for different regions in the Dalitz plots. The cusp structure at 
the NE threshold is described with the Flatte formalism and its variation in the Dalitz plot is analyzed. 


1 Introduction 

The associated strangeness production of kaon and lambda 
is the energetic lowest possibility of creating particles with 
open strangeness in nucleon nucleon collisions. This ele¬ 
mentary reaction is of fundamental interest as it involves 
the dissociation of an ss quark pair into two hadrons with 
s = 1 and s = —1. 

The creation or dissociation of the ss quark pair is de¬ 
scribed in different models by the exchange of a meson 
between the incoming nucleons. In most models kaon or 
pion exchange is assumed, but this approach can be ex¬ 
tended to 77, p, and oj mesons and to the strange K* mesons 
. To separate the contributions of the strange 
and non strange meson exchange experimentally, two pos- 
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sibilities are proposed: First significant differences of the 
analyzing power and the spin transfer coefficient for both 
exchanges are predicted [5]. Secondly the evidence of nu¬ 
cleon resonances in the lambda kaon subsystem can only 
be explained by non strange meson exchange [6]. 

The associated strangeness production in proton pro¬ 
ton collisions has been extensively studied by the DISTO 
collaboration [7], the GOSY-TOF collaboration [5], [5], 
close to threshold by the GOSYll collaboration Bn], and 
at a higher beam momentum by the HADES collabora¬ 
tion m- In spite of this large experimental data base 
no conclusive solution has been found on the exchange 
mechanism. While measurements of the spin transfer co¬ 
efficient of the DISTO experiment indicate kaon exchange, 
COSY-TOF measurements yield a large N* contribution 
[8], which is a sign for dominant pion exchange. 

Apart from the reaction mechanism itself, the associ¬ 
ated strangeness production offers powerful tools to study 
a) the decay of N* resonances into channels with strange¬ 
ness, b) the hyperon nucleon interaction, and c) the cou¬ 
pled channel effect of (NA NE) at the NE threshold. 

a) The first evidence for the creation of N* resonances 
in the associated strangeness production were found by 
comparing the Dalitz plot distributions measured with the 
COSY-TOF experiment [H1B2 with predictions, which as- 
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sume a strong impact of the resonances N(1650), N(1710), 
and N(1720). Further insight to the weight of these reso¬ 
nances and their branching ratio to the KA channel can 
be obtained by a partial wave analysis, with precise data - 
including polarization observables - as input. Therefore, 
a combined analysis of all existing data with the Bonn- 
Gatchina partial wave analysis program m has been star¬ 
ted [13] . First results of the partial wave analysis of HADES 
data are given in m- 

b) The invariant mass distributions of the associated 
strangeness production allows the interaction of the in¬ 
volved particles to be studied as a function of the relative 
energy in the range from zero to several hundred MeV. Es¬ 
pecially the measurement of the nucleon hyperon invariant 
mass distributions starting from zero interaction energy 
enables the determination of the nucleon hyperon scatter¬ 
ing length [Bin]- The effective pA scattering length has 
been determined from COSY-TOF data [IS]- A precise 
knowledge of the N* resonances is needed, as they influ¬ 
ence the pA invariant mass spectrum by reflections of the 
KA channel, thus limiting the precision of the scattering 
length determination. 

c) The last tool provided by the associated strangeness 
production is to study the coupled channel effect at the 
threshold of NE in the nucleon hyperon invariant mass 
distribution. The cusp structure at the threshold was dis¬ 
covered in 1961 by a bubble chamber measurement of K“ 
absorption in deuteron leading to pAtt" [T9]- This was 
confirmed by higher statistic measurements of the same 
reaction mm- The measurements were described by dif¬ 
ferent models and calculations [221I231I24] . As for the pKA 
final state both thresholds of NS and KE have to be ex¬ 
ceeded, the cusp effect can be observed for beam momenta 
above 2.7GeV/c. First evidences for a cusp in this final 
state at beam momenta of 2.75GeV/c and 2.85GeV/c 
were reported of measurements with the GOSY-TOF ex¬ 
periment m- A detailed examination of the cusp was 
performed with data measured at a beam momentum of 
3.05GeV/c [25]. 

In addition to the fundamental processes which can be 
studied with the associated strangeness production, this 
elementary reaction is the basis of more complex reac¬ 
tions, which are under investigation and which need the 
precise knowledge of its observables. Among these com¬ 
plex reactions are the properties of hypernuclei [B; tti® 
possible existence of ppK cluster or bound states EH, and 
the prediction of strangeness enhancement in heavy ion in¬ 
teractions as a sign for the quark gluon plasma [55] . Even 
information of the associated strangeness production mea¬ 
sured at other kinematic regions can contribute to a better 
understanding, as they can help to describe the involved 
nucleon resonances and threshold effects, which have to be 
subtracted in order to identify these higher order effects. 

The purpose of this paper is to present the analysis of 
a high precision measurement of the reaction pp ^ pK+A 
in order to increase the existing data base in this momen¬ 
tum region by a factor of 5. Special emphasis is placed 
on the description of the threshold effect. Angular distri¬ 
butions, Dalitz plots, and invariant mass distributions are 


given and examined. The polarization variables will be 
presented in a forthcoming paper. 

2 Experimental setup 

The GOSY-TOF experiment is a non-magnetic spectro¬ 
meter, which is situated at an external beam-line of the 
accelerator GOSY in the research center Jiilich. It consists 
of a cylindrical vacuum vessel of 3.5 m length and 3 m dia¬ 
meter (see hg. [T]) . The pressure inside the vessel is lower 
than 6 x 10“^ hPa. 



Fig. 1. Side view of the COSY-TOF spectrometer. In beam di¬ 
rection the start counter (Start), the straw tube tracker (STT), 
the barrel scintillators, the inner ring (Quirl), the outer ring 
(Ring), and the Calorimeter are shown. All detectors and the 
liquid hydrogen target are located inside the vacuum vessel. 


2.1 Scintillation trigger detectors 

At the inner walls of the cylinder part of the vessel 96 plas¬ 
tic scintillator bars are attached, which have a thickness 
of 15 mm, a width of ca. 10 cm and a length of 285 cm. 
They are arranged to form a slightly conical barrel of 
about 300 cm in diameter. The barrel scintillator is de¬ 
scribed in ref. [53]. The downstream end-cap of the ves¬ 
sel is covered by three layers of 5 mm thick plastic scin¬ 
tillators, which are arranged in two rings, the inner one 
ranging from 2.4 cm to 58 cm radius, the outer one from 
54 cm to 150 cm. One of the layers is composed of 48 (in¬ 
ner ring) and 96 (outer ring) wedge-shaped scintillators. 
Each of the other two layers consists of 24 (inner ring) 
and 48 (outer ring) elements, formed as Archimedian spi¬ 
rals. The spirals change the direction of rotation from one 
to the other layer, therefore, 1104 (inner ring) and 2304 
(outer ring) pixels can be defined by coincidences. The in¬ 
ner ring is described in ref. [30] . The scintillators are read 
out by XP2020 photo-multipliers. Downstream, the inner 
ring counter is followed by a scintillator calorimeter. The 
data of the calorimeter are not used for this measurement. 
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2.2 Target and beam defining detectors 

The liquid hydrogen target is located in the entrance part 
of the vacuum vessel. It is a cylindrical cell with a dia¬ 
meter of 6 mm and a length of 4 mm, which consists of 
ss 4 pm thick copper walls in the barrel part and of 0.9 pm 
thick Mylar window foils in the entrance and exit. In or¬ 
der to reduce rest gas deposition onto the target windows, 
the pressure surrounding the target cell is kept below 
2 X 10“^ hPa. A scintillation start counter is located at a 
distance of about 15 mm downstream to the target center. 
It is made from two 1 mm thick discs, with an outer ra¬ 
dius of 75 mm and a central beam hole of 1.5 mm radius. 
Each disk is divided into 12 wedge-shaped parts, which 
are individually read out by Hamamatsu R1450 photo¬ 
multipliers. A double-sided silicon detector disk with a 
thickness of 300 pm and a central beam hole of 2.75 mm 
radius is located 26 mm downstream. The data of this de¬ 
tector are only applied for pKE+ analysis, it will be de¬ 
scribed in a forthcoming paper on the E"*" results. In order 
to veto beam particles which hit directly the start counter, 
a scintillation counter with a hole of 1 mm radius is in¬ 
stalled 50 mm upstream of the target. Further upstream 
at 50 cm and 100 cm two veto counters with central holes 
of 6 and 15 mm diameter detect beam halo particles which 
are around the beam core. The total beam rate and its 
X- and y-intensity distribution is measured with a fiber 
hodoscope located 1 m downstream of the spectrometer 
exit. The target, start and veto counters are aligned to 
the center of the beam line with a precision in the order 
of 0.1mm with an optical telescope, which can be placed 
2 m upstream of the target into the beam line. The beam 
is focused onto the target with a profile which has a full 
width at half maximum of about 1 mm. Due to precise 
tuning of the beam line magnets the maximum rate in the 
veto counters is less than 1% of the total beam rate. 


2.3 Tracking detector 

The main detector system is the central straw tube tracker 
(STT). The STT is placed about 24 cm behind the target 
inside the vacuum tank. It consists of 2704 straw tubes, 
which are arranged in 13 double layers and fixed in six 
orientations with an angle of 60° to each other in order 
to enable 3D track reconstruction. A single straw has a 
length of 1050 mm, 10 mm diameter and 30 pm wall thick¬ 
ness. The active detector volume of the 13 double lay¬ 
ers is about 1 m in diameter and 24 cm in depth. An in¬ 
ner beam hole in each double layer with 15 mm diameter 
avoids beam interactions with the detector. Each straw 
tube is made of aluminized Mylar film and filled with Ar 
-I- CO 2 (80 : 20) gas mixture at a pressure of 1.2 bar. 
The straws are self-supporting by the over pressure of the 
chamber gas. Therefore, sufficient mechanical stability is 
provided, despite the low material budget of A/Aq « 1% 
for 26 layers of straws. The anode is a 20 pm thick gold 
plated tungsten-rhenium wire which is stretched along the 
straw axis and held at a potential of -1-1820 V m- 


The arrangement of the target and detectors facilitates 
a volume in the z-direction from the target to the begin¬ 
ning of the straw detector, which is - apart from the start 
counter and silicon counter - free of any material. There¬ 
fore, a secondary vertex in this volume indicates with high¬ 
est probability a A or a Kg decay. 

2.4 Trigger 

The trigger is defined by signals of the start counter, the 
barrel scintillators and the end-cap scintillators. At least 
one signal of the start counter and at least 4 hits in the 
barrel and end-cap are required. A hit in the ring counters 
is defined by coincident signals in two of the three layers 
and a hit in the barrel counter is defined by a coincidence 
signal of the photo-multipliers, which are situated on both 
ends of a barrel scintillator bar. The coincidence width is 
100 ns, which takes into account the different time of flight 
of particles between the start and stop counters and of 
different light paths in the scintillators and light-guides. 
In order to facilitate the rate of « 2000 events/s written 
to disk, only the STT data were read out, as this read 
out electronics has a dead time of 10 ps. The dead time of 
the ADC and TDC readout for the scintillators is in the 
order of 200 ps. Therefore, data which include the ADC 
and TDC information of the scintillators are recorded for 
controlling the trigger system for only about 5% of the 
beam time. 


3 Analysis 

3.1 Event Reconstruction 

The analysis of the pKA events is based on the tracking 
information obtained from the straw tube tracker. The de¬ 
termination of the correlation between drift time and the 
track distance to the anode is performed in two steps [32l 
133] : In the first step the correlation is obtained by an in¬ 
tegration over all measured drift times and by assuming 
a homogeneous illumination of the straw tube by particle 
tracks. In the second step the position in the straw tube 
is calculated from reconstructed tracks, which are cross¬ 
ing this tube. With this information the correlation be¬ 
tween drift time and distance to the center is recalculated. 
The spatial resolution of the STT is cr = (140 ± 10) pm at 

2.5 mm distance to the anode |33) . 

The main part of the program which extracts the tracks 
from the straw tube isochrones is contained in the tof- 
Straw library [M] and described in detail in [33] , [35] , and 
m- The primary and decay vertices are submitted to a 
geometrical fit, which includes the information, that the A 
decay plane has to contain the primary vertex. The events, 
for which the fit converges, are submitted to a kinematic 
fit by taking the geometrical fit results as start values. The 
kinematic fit applies momentum and energy conservation, 
the masses of the A and its decay particles are directly 
used as input, as the A is identified by the secondary ver¬ 
tex and the decay proton has always the smaller angle to 
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the A direction. The assignment of the proton and kaon 
masses to the primary tracks is done by applying the kine¬ 
matic fit for each of the two solutions, the result with the 
lower is chosen. Events with a x^/NDF of less than 5, 
a decay length in the rest system of the A of larger than 
2 cm and an angle of the decay proton to the A direction 
of larger than 2.5° are evaluated as pKA events. 

3.2 Monte Carlo calculations 

The corrections for the detector acceptance and the recon¬ 
struction efficiency are done with a Monte Carlo program 
(GEANT 3.2 [37]) with phase-space distributed events as 
input. The corrections were controlled with a Monte Carlo 
program with an input, which simulates the measured 
proton cm angular distribution. No significant changes of 
the resulting angular distributions were found. The com¬ 
bined acceptance and reconstruction efficiency of pp —>■ 
pK(A ^ p7r“) amounts to (22.5 ± 1) %. The primary ver¬ 
tex resolution is (7(x, y) = (0.5 ±0.1) mm and a{z) = 
(1.5 ± 0.2) mm. The transverse resolution of the secondary 
vertex is in the same range as of the primary one, while the 
longitudinal resolution is worse: a(z) = (2.0 ± 0.5) mm. 
The resolution of the pA invariant mass is (j(mpA) = 
(1.1 ± 0.3)MeV/c^, the resolution of the KA and pK in¬ 
variant masses are cr(mKA) = (1-1 + 0.3)MeV/c^ and 
o'(wpK) = (1.0 ± 0.3) MeV/c^ [SSIISB] . 

3.3 Background 

One source of background are multi pion events. They can 
produce spurious secondary vertices if tracks are recon¬ 
structed erroneously. These spurious secondary vertices 
are expected to be close to the target, as the density of 
tracks is highest in this region. These background events 
can be identified as deviation of the measured from the 
calculated A decay length (fig. |2|). The normalized data 
are htted with /(x) = po ■ exp(— x/ctq) in the range be¬ 
tween 4 cm and 30 cm. po is the ht coefficient and ctq is the 
literature value of the A decay length (7.89 cm) [35]. The 
background events are concentrated in the first centime¬ 
ters of the decay length. The excess of the data compared 
to the prolongation of the fit curve between 2 cm and 4 cm 
adds up to 5% of the total events. 

The background events caused by > yA cannot 
be detected by geometry information, as the event topol¬ 
ogy is similar to the pKA events. The suppression of these 
events by the cut of the kinematic fit is studied with 
a Monte Carlo sample of pp ^ pKE°(^ Ay). The re¬ 
construction efficiency of this sample is 8.1%, which has 
to be compared to the reconstruction efficiency of pp ^ 
pKA which amounts 22.5%. Including the ratio of the 
cross sections of 3.1pb/21.8pb a contamination of 5% 
of events is obtained. The analysis of the Monte Carlo 
data of pKS° shows, that these events are shifted by the 
kinematic fit to the backward region of the A and kaon 
cm scattering angles. For pKA events this effect is not 
observed. 



Fig. 2. A decay length distribution in the A rest system. The 
data are corrected for acceptance and efficiency. The fit with 
Po ■ exp(—a;/7.89 cm) is applied in the range between 4 cm and 
30 cm (red line). The prolongation between 2 cm and 4 cm is 
shown with the blue line. 


4 Results 


4.1 Angular distributions in the cm system 


The angular distributions in the center of mass system are 
shown in fig. |3| The data are corrected for the detector ac¬ 
ceptance and reconstruction efficiency and are normalized 
to the total cross section of 21.8 pb, which is determined 
in a former COSY-TOF measurement |^. As target and 
projectile are identical particles, which cannot be distin¬ 
guished in the cm system, these distributions have to be 
symmetric to cos'i?°“ = 0. Small deviations from symme¬ 
try are existing for all three distributions: At backward 
angles there is a surplus of events, which is mainly due 
to the pK(S° ^ Ay) background. At extreme forward 
angles (cosi?°™ > 0.95) of the proton and kaon distribu¬ 
tions a deficit of events is found. The affected range cor¬ 
responds to scattering angles in the lab system of smaller 
than 5° - 6° degrees. It is supposed that the deficit of 
events is caused by an efficiency drop in the cones of the 
start counter, which cannot be exactly described by the 
Monte Carlo program. 

While the angular distribution of the heavier reaction 
products, p and A, exhibit distinct deviations from uni¬ 
formity, the distribution of the lighter kaon is nearly con¬ 
stant. This behavior is examined in detail in section EH 
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Fig. 3. The angular distributions in the cm system are shown 
for the proton (upper figure), kaon (middle), and A (bottom). 
The results from a former COSY-TOF measurement [S] are 
given as red squares. Beneath each distribution the detector 
acceptance and reconstruction efficiency e is shown. 


Fig. 4. The angular distribution in the pA Gottfried-Jackson 
frame is shown in the upper figure. Middle: pK Gottfried- 
Jackson frame, bottom: KA Gottfried-Jackson frame. The re¬ 
sults from a former COSY-TOF measurement are given as 
red squares. Beneath each distribution the detector acceptance 
and reconstruction efficiency e is shown. 


4.2 Gottfried-Jackson angular distributions 

The angular distributions in the three two particle subsys¬ 
tems of the final state pKA are constructed by boosting 
the momenta into the rest system of the two particle sys¬ 
tem. The angle between the boosted beam vector and one 
of the particles of the two particle system is the Gottfried- 
Jackson (GJ) angle. It is denoted with where i, j are 
the particles of the two particle rest system and b stands 
for the boosted beam direction. The GJ angular distribu¬ 
tions are shown for the rest systems of Ap, pK, and KA in 


fig. m In contrast to the cm distributions the GJ angular 
distributions do not have to be symmetric to cos'dJ;J; = 0. 

While only in the KA rest system higher angular mo¬ 
menta due to nucleon resonances are expected, the GJ 
angular distributions exhibit the highest anisotropy for 
the Ap rest system. The GJ angular distributions are ex¬ 
amined in dependence on bins of the invariant masses in 
section 15.2.21 
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4.3 Dalitz plots 

The Dalitz plots are normalized by scaling the maximum 
channel content to 1. The Dalitz plot of versus rn^j^ 
is shown in fig. [S] Two dominant structures, which are 
stretched mainly in the vertical direction, are observed. 
The first one is located at low pA invariant masses and 
arises from the pA final state interaction [3]. The second 
structure located at the NE threshold has its maximum in¬ 
tensity around TOka 2.85GeV^/c'^ and extenuates with 
rising KA invariant mass. 

The Dalitz plot of versus rn^^ is shown in fig. [51 
The two dominant structures of the previous Dalitz plot 
are emerging as diagonal spread elevations. Again the struc¬ 
ture connected to the NE threshold exhibits a strong in¬ 
tensity variation across the Dalitz plot. 

No structures at the KE threshold or near the reso¬ 
nance masses can be detected. A detailed analysis of the 
invariant mass distributions is given in section 05.21 


5 Detailed examinations and discussions 

5.1 CM angular distributions 

Former COSY-TOF measurements exhibit that the / = 2 
Legendre polynomial coefficient, which is a measure for 
the contribution of higher partial waves, rises with increas¬ 
ing beam momentum The data presented here allow 
for the first time to extract the Legendre coefficients not 
only from the integrated angular distribution but from 
a series of angular distributions, each constructed for a 
bounded range of the center of mass momentum. 

For all reaction particles (p, K, A) the magnitude of the 
center of mass momentum of the particle is divided into 

6 equidistant intervals. For each momentum interval the 
angular distribution of the particle has been fit according 
to the formula 


-5 -= C ■ (apPo + a 2 P 2 (cosi?™) a 4 P 4 (cosi?™)) 

d cos 

.( 1 ) 

C is a normalization constant. This formula describes 
the differential cross section directly as a sum of Legen¬ 
dre polynomials Pi by neglecting the interference terms. 
Due to identical particles in the initial state, the angular 
distributions in the center of mass system have to be sym¬ 
metric with respect to cos'd'^™ = 0. Therefore, the coeffi¬ 
cients oi and 03 are set to zero. With this parametrization 
and the restriction to the maximum angular momentum of 
I = 2, the coefficient 02 contains products of SD, PP, and 
DD partial wave amplitudes, and 04 contains only the DD 
partial wave product. The angular distributions are nor¬ 
malized for each range of the center of mass momentum 
separately, in order to be independent of the variation of 
the cross section with the center of mass momentum. 

The angular distributions are shown in the appendix 
(fig. [13] for the proton, fig. [20] for the kaon and fig. [2T] for 
the A particle). In addition in each figure the Legendre 
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Fig. 5. Top: The Dalitz plot of the measured data. In order 
to enhance the visibility of the structures the color code is 
truncated to 0.85 (less than 0.1% of the bins have an occupancy 
of larger than 0.85). The black horizontal line indicates the KE 
threshold, the vertical one the NE threshold. The red area at 
the ordinate indicates the mass of the N*(1720) resonance, the 
green area indicates the N*(1710) resonance, and the blue area 
the N*(1650) resonance. Middle: The projection on the squared 
invariant mass pA is shown. This projection is corrected by 
the combined detector and reconstruction efficiency, which is 
shown in the bottom frame. The NE threshold is indicated by 
the black line. 


polynomials Pp, P 2 , and P 4 , weighted with the correspond¬ 
ing fit coefficients, are plotted and the corresponding mo¬ 
mentum range is indicated. For the proton and kaon the fit 
is restricted to the range of —0.95 ^ cosd™ ^ 0.95. The 
limit in the backward direction takes into account the in¬ 
crease of background at this region. The cut in the forward 
range corresponds to a limit in the laboratory scattering 
angle of about 5°, and takes into account that the pro¬ 
ton and kaon have to be detected by the start-counter, 
of which the efficiency determination for the very small 
angles has a larger error compared to other angles. 

The variations of the coefficients ap, 02 , and 04 with 
the center of mass momentum are shown in figs. [71111111 for 
p, K, and A, respectively. The coefficient 02 of the proton 
rises nearly linearly from 0.15GeV/c to 0.35GeV/c and 
is constant in the momentum region 0.35 — 0.5GeV/c. In 
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Fig. 6. Top: The Dalitz plot of the measured data. In order to 
enhance the visibility of the structures the color code is trun¬ 
cated to 0.85 (less than 0.1% of the bins have an occupancy 
of larger than 0.85). The black vertical line indicates the KE 
threshold. The red area at the abscissa indicates the mass of 
the N*(1720) resonance, the green area indicates the N*(1710) 
resonance, and the blue area the N*(1650) resonance. Middle: 
The projection on the squared invariant mass KA is shown. 
This projection is corrected by the combined detector and re¬ 
construction efficiency, which is shown in the bottom frame. 
The KE threshold is indicated by the black line. 


contrast to this behavior the coefficient 02 of the A rises 
with the third power of the cm momentum: 

02 = 1/C • / • (pr[GeV/c])3 

with / = (12.2 + 1.4) ^b/(GeV/c)^ (see red line in 
fig. ini). C is the normalization factor defined in equation 

m 

The coefficient 04 of the proton starts to be signifi¬ 
cantly different from 0 at the cm momentum of 0.3 GeV/c 
and rises up to 1/C • 0.6 [fib/sr] at 0.45 — 0.5 GeV/c. The 
fit of the A distribution gives only for the momentum bin 
of 0.5GeV/c a significant 04 contribution with about 1/3 
compared to the proton angular distribution. The kaon 
angular distribution can be described for all momentum 
bins by the coefficient oq alone, no higher Legendre poly¬ 
nomials are needed for the fit. The reduced of the fit 
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Fig. 7 . The dependence of the fit coefficients of the proton 
angular distribution on the proton center of mass momentum 
is shown. The upper small figure shows the corresponding re¬ 
duced of the fit. 



kaon cm momentum [GeV/c] 


Fig. 8. The dependence of the fit coefficients of the kaon an¬ 
gular distribution on the kaon center of mass momentum is 
shown. The last bin is missing, as the maximum kaon momen¬ 
tum is less than 0.44GeV/c. The upper small figure shows the 
corresponding reduced of tho fit- 
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A cm momentum [GeV/c] 


Fig. 9. The dependence of the fit coefficients of the A anguiar 
distribution on the A center of mass momentum is shown. The 
upper small figure shows the corresponding reduced of the 
fit. The distribution of the coefficient 02 is fitted with c • 
(continuous red fine). 


increases only by a factor of 1.1 to 1.3 by omitting the P 2 
and P 4 polynomials in the fitting procedure. 

The isotropy of the kaon angular distribution, even for 
the highest momenta, is in contrast to the behavior of 
the uj in the pp —> ppw reaction [39]. Here the cm angu¬ 
lar distribution of the lighter particle (w) is strongly an¬ 
isotropic (for the measurement with a beam momentum of 
3.059 GeV/c the maximum cm momentum of the to is with 
0.4GeV/c approximately in the same range of the maxi¬ 
mum cm momentum of the kaon). The protons exhibit a 
nearly flat angular distribution in contrast to the proton 
or A of the pKA reaction. A better comparison is the pp 
pp ?7 reaction, as the mass of the 77 , compared to the one of 
the oj, is closer to the kaon mass. Data measured at a beam 
momentum of 2.94GeV/c, where the rj has 0.554 GeV/c 
cm momentum, yield a strongly anisotropic angular dis¬ 
tribution. In contrast, the angular distribution is isotropic 
at a beam momentum of 3.67GeV/c, corresponding to 
7] cm momentum of 0.774 GeV/c [40]. Therefore, there is 
no clear sign, whether the cm angular distribution of the 
lightest particle of a three particle system is particularly 
isotropic due to kinematical effects. 

5.2 Invariant mass distributions 

For better comparability with literature - especially con¬ 
cerning the threshold examinations - the linear invariant 
mass distributions are analyzed. 

The measured invariant mass distribution is given by 


= PS(mij) ■ |M(mij)p • e(mij) (2) 

PS denotes the phase-space volume, M(mij) the re¬ 
action amplitude, and e is the detector acceptance and 
reconstruction efficiency, which are dependent on the in¬ 
variant mass. By modeling the reaction with the Monte 
Garlo technique the phase-space, the acceptance and re¬ 
construction efficiency are calculated: 

du^C 

= PS(mij) •e(mij) (3) 

The division of the measured distribution by the Monte 
Garlo distribution generates a spectrum, which is only de¬ 
pendent on the reaction amplitude M(mij) and no assump¬ 
tions on the normalization to the phase-space distribution 
are neededQ In the following these spectra are analyzed. 

5.2.1 Proton-A invariant mass 

The main goal of the examination of this invariant mass 
distribution is the analysis of the behavior at the NE 
threshold. In the strict sense there are two thresholds: nE+ 
at 2.I29GeV/c^ and pE° at 2.I3IGeV/c^. This mass dif¬ 
ference is about Icr of the invariant mass resolution and 
not resolved in the data. In the following calculations the 
mass of the lower threshold is used, as this channel dom¬ 
inates the cusp effect [35]. For convenience the threshold 
is named NE threshold. 

In order to examine the NE threshold we assume as a 
working hypothesis that the final state interaction (FSI), 
the NE threshold (TH), and the reflections of N* reso¬ 
nances (RF) are independent and interference terms are 
negligible: 

do-meas 

—-/-- = FSI(topa) + TH(mpA) + RF(mpA) (4) 

dmpA drUpA 

Of course this is a simplification, as the reflections of 
the N* resonances can interfere with the final state interac¬ 
tion and the threshold effect. The final state interaction, 
the reflection of N* resonances and the threshold effect 
are described with functions with free parameters. The 
parameters are determined by fitting the invariant mass 
distribution with the sum of these functions. 

The final state interaction is treated like the other 
terms as an additive term in order to incorporate that 
the FSI has strong influence at low pA invariant mass 
and this influence is negligible towards larger invariant 
masses0 This is in contrast to the ansatz, where the reac¬ 
tion amplitude is factorized with the FSI amplitude (see 
for instance [3] and references therein). 

^ In publications of the DISTO Collaboration (e.g. | 41| ) this 
distribution is named “deviation spectrum”. 

^ A similar behavior is measured for the pp final state in 
pp —* ppp at an excess energy of 41 MeV |43| and in np —> ppn~ 
at excess energies between 405 and 465 MeV [43] . 
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The final state interaction is parameterized with 


FSI(mpA) = exp(co H--) (5) 

WpA - C2 

The effective scattering length can be expressed by two 
of the three coefficients of equation [5J 


Oeff = —1/2 hcciA 




TOpTOA - cs) • (m^ 


C 2 )" 


( 6 ) 

Winin and TTimax denote the range of the invariant mass 
nipA in which this formula is applied to determine the scat¬ 
tering length. Equations [ 5 ] and [6] are given in [35], they are 
derived in [16] . Equation [5] is applied to express the co¬ 
efficient Cl with the effective scattering length. The value 
of tteff = 1.233 fm is taken, which is deduced from data 
measured at a beam momentum of 2.7 GeV/c [35]. At this 
beam momentum no disturbing effects of the N27 thresh¬ 
old are present. The influence of the N* resonances were 
tested by evaluating different regions of the Dalitz plot, 
they give limits of the measured effective scattering length 
between -1.4 and -1.2 fm [35] . at the beam momentum of 
2.95 these limits were found to be larger with -0.9 and -2.1 
fm [18]. With this procedure the final state interaction at 
2.95GeV/c is expressed with two parameters cq and C 2 . 

The nucleon resonances N*(1650), N*(1710), N*(1720) 
contribute to deviations from the phase-space in the Dalitz 
plot [5]. As these resonances decay into K+A their in¬ 
fluence can be seen directly in the invariant mass distri¬ 
bution of tok+a as reflections in the invariant mass 
distribution of rUpA- As their widths are in the range of 
50 — 400 MeV [35] the reflections are slowly varying with 
TTipA compared to the ESI and the threshold effect, whose 
structures are smaller than 20MeV/c^. Therefore, these 
reflections are parameterized as a polynomial of second or¬ 
der. The main purpose of this parametrization is to quan¬ 
tify the background beneath the cusp effect, and not to 
draw conclusions on these resonances. For simplicity this 
part is named “reflections” in the further description. 

As the threshold effect is expected to be a coupled 
channel effect, we apply the Flatte formalism [4^. This 
formalism was developed to describe the rjir (tttt) inter¬ 
action at the two kaon threshold. It is assumed that the 
ao(980) (/o(980)) resonance changes its width due to the 
opening channel of KK thus inducing a cusp structure. 
This formalism is adapted to the pA channel by assuming 
a virtual state of pA, which decays into pA and into NE, 
thus causing a cusp effect E3- 


TH(mpA) = C\ 


^ta/^0 * Fap 


to ; 


^ ^ im^{rAp + Ten) 


(7) 


TOi- and To are the mass and the total width of a virtual 
pA state. C is a normalization parameter. As it cannot be 
separated from the unknown total width Tq, the product 
C-Tq is treated as one parameter for the fitting procedure. 
Tap and Ten are the partial decay widths of the virtual 
state, which are given by: 


Tap = ffAp • Q'Ap and Ten = ffSN ■ 9 en (8) 

gAp and (;en are the coupling constants of the two 
channels to the virtual pA state and qAp and qen are the 
cm momenta in the corresponding two-body subsystems. 
In dependence on the invariant mass TOpA, (Jen is imagi¬ 
nary below the NE threshold: 


9en = I 


^ ((toe+ + TO„)2 - TO^^) • (to2^ - (toe+ - TOn)2) 


and real above the NE threshold: 


(9) 


(Jen = 


\l ("IpA - (™s+ + • (WpA - (r«E+ - 


2TOpA 


( 10 ) 

The cm momentum qAp is always real and given by ex¬ 
changing the E”*" and neutron mass (toe+ , TOh) with the A 
and proton mass (toa, rup) in equation [TUI A detailed ex¬ 
amination of Flatte-like distributions [4^ reveals that the 
coupling constants and the mass of the virtual state show 
a scaling behavior. Therefore, the coupling constant gAp 
is fixed to the arbitrarily chosen value of 1, and only the 
remaining parameters CTo, ^en, and TOr are determined 
together with the parameters for the reflections and FSI 
by a fit procedure. 

The measured invariant mass distribution - divided 
by the Monte Carlo distribution - fitted with the function 
as described in equation |3| is shown in fig. [TU] The fit 
starts at the kinematic limit of rUp + toa and is limited 
to 2.250 GeV/c^, which is 8MeV/c^ below the maximum 
possible kinematical value. This limit is chosen as the data 
and Monte Carlo distributions have above this value a 
sharp decrease to zero, therefore, the uncertainties of the 
ratios are large. 

The structures of the data are well described by the fit. 
The sharp rising leading edge, the peak at the threshold, 
and the slow decreasing above the threshold is reproduced 
by the Flatte formalism. The mass of the virtual pA state 
is 2.31GeV/c^. However, this value is depending on the 
arbitrarily chosen value of the pA coupling constant to 
this virtual state. Changing this coupling constant for in¬ 
stance from 1 to 0.5 produces a fit result for the mass of 
the virtual pA state of 2.20GeV/c^. As remarked in |45j . 
the mass of the virtual state has to be larger than the 
threshold in order to generate a sharp cusp shape. For 
masses below the threshold, the distribution varies more 
smoothly, and is more similar to a Breit-Wigner distribu¬ 
tion. The coupling constant of the virtual state to the NE 
channel ((Jen = 1.253 ± 0.189) is larger than the coupling 
constant to the pA ((jAp = Ij fixed value) channel. This 
is in contrast to the previous analysis |25| . which applied 
a different mass of the virtual state. The first attempt to 
describe the threshold behavior of the pA invariant mass 
distribution with the Flatte formula was performed in 
where a measurement of K d ^ pAtt is described. Here 
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Fig. 10. The whole data set divided by the Monte Carlo data 
is shown as dots with error bars. The fit result is given by the 
red line, it is composed of the FSI part (blue line), the part of 
the reflections (black line), and the Flatte-distribution (green 
line). The NE threshold is marked with the vertical black line. 
The kinematic limits are indicated by two dashed lines. 


the best fit result is obtained with a mass of 2.20GeV/c^ 
as input parameter. 

As can be seen in the Dalitz plot of fig. [SI the density 
in the NS threshold region, is strongly decreasing with in¬ 
creasing rriKA- In order to examine, whether this variation 
is due to the cusp effect or due to the N* contributions, 
the Dalitz plot is divided into 6 slices, by applying cuts 
on the helicity angle between the proton and the kaon in 
the pA subsystem (s. fig. HH). The limits of each slice are 
given by 


- 1 -h z • 1/3 < coszJp^"^ < -1 -h (z -f 1) • 1/3 (11) 

with i = 0..5. The projection of these slices onto the 
pA invariant mass results in identical distributions if the 
Dalitz plot is filled according to phase-space. 

For each slice the invariant mass pA is fitted with the 
same fnnction as applied for the whole data set. As the 
FSI distribution is expected to be independent of the he¬ 
licity angle, the parameters cq, C2 are fixed to the values 
which are obtained by the fit of the whole data set. In 
addition it is assumed that the cusp does not change its 
shape, therefore, the parameters mj-, and are fixed 
to the values obtained by the whole data fit. CTq and 
the coefficients of the reflections are determined for each 
slice with the same fit procedure. The height of the cusp 
is proportional to CTo. With nir = 2.31GeV/c^ the pro¬ 
portional factor is 1.39. As the shape of the cusp is con¬ 
stant, CTo is in addition proportional to the area beneath 



Fig. 11. The division of the Dalitz plot into six slices by se¬ 
lecting ranges of the helicity angle between proton and kaon in 
the pA subsystem is shown. If the Dalitz plot is filled according 
to the phase-space, the projection of the slices i = 0..5 onto 
the invariant mass pA results in identical distributions. 


the cusp described by the Flatte function. The amount 
of the reflections beneath the cusp is determined as the 
value of the second order polynomial at the cusp posi¬ 
tion (zTZpA = 2.129 GeV/c^). The individual distributions 
of the slices and their fit results are shown in the appendix 
(fig. m]). The variation of the height of the Flatte distri¬ 
bution and the amount of the reflections at the threshold 
is shown in fig. [121 While the height of the Flatte distribu¬ 
tion varies only slowly with the helicity angle, the part of 
the reflections is nearly constant for the first three slices 
and then drops from 0.3 to 0 in the last three slices. This 
indicates that the large intensity variation in the Dalitz 
plot along the cusp line is caused by the background gen¬ 
erated by the N* resonances. 

The assumption, that the cusp does not change the 
shape, seems to be justified by the ht results. It can be 
seen in fig. [521 that the cusp is well described by the fit 
function, having the same parameters, which affect the 
cusp shape. The quality of the fit is nearly the same for 
each slice, all reduced are in the range of 1.4 ± 0.2. 

The angular distribution of the cusp effect reflects the 
angular distribution of the pA system, which stems from 
the NE system at the threshold. The angular momenta 
between p and A inside the pA system are described by 
the Gottfried-Jackson angular distribution (coszJp^^). As 
the pA system has the same spin-parity as the NA system 
at the threshold (J^ = 0+, I"*"), it can only be in relative 
S or D waves |25j . The angular momentum of the cusp can 
be examined with the kaon angular distribution in the cm 
system, which is one of the two particle subsystem K-(NZ’) 
and thus the kaon angular distribution is the mirrored one 














The COSY-TOF Collaboration: High statistics measurement of the reaction pp pK+A 


11 




2 - 

Y 


cusp 
N* reflections 



Fig. 12. The height of the cusp (green crosses) and the 
amount of the reflections at the cusp position (black crosses) 
for the different regions of the Dalitz plot as shown in fig. 1111 
The blue dashed line shows the fit result of the height of the 
cusp of the total data set. As the errors of the fit parameter 
CTo are found to be asymmetric, the errors of the MIGRAD 
routine of the ROOT | 47| program are shown. The upper small 
figure shows the corresponding reduced of the fit. 

of (NA'). This angular distribution has to be symmetric 
to cosYg" = 0. Both angular distributions cosY™ and 
cosYp^^ of the height of the Flatte distribution and the 
amount of the reflections beneath the cusp are generated 
with the slice method as described in the previous section. 
In this case cosY™ and cosYp^^ are divided into 6 bins. 
The fit results for the cosY™ slices are shown in fig. [231 
and for the cos Yp^^ slice in fig. [24] in the appendix. 

The variation of the height of the Flatte distribution 
and the amount of the reflections at the threshold in de¬ 
pendence on the different slices is shown in figs. [T31 dH 

No systematic deviations of the cusp height distribu¬ 
tions from isotropy, which could be a hint for P or D waves, 
can be seen. This is a confirmation of the findings in [23] , 
which are based on a different data set of COSY-TOF 
and on a different analysis method. The shape of the cusp 
does not change with the examined angles, as can be seen 
from the fits of the individual slices. The reduced of all 
fits does not deviate from the mean value of 1.36 by more 
than 0.2. 


5.2.2 Kaon-A invariant mass 

Deviations of the KA invariant mass from phase-space can 
be caused by a KE threshold effect or by nucleon res¬ 
onances, which decay into the KA channel. In addition 


Fig. 13. The height of the cusp (green crosses) and the 
amount of the reflections at the cusp position (black crosses) 
for the different regions of cos Y™. The blue dashed line shows 
the fit result of the height of the cusp of the total data set. As 
the errors of the ht parameter CTo are found to be asymmetric, 
the errors of the MIGRAD routine of the ROOT |47| program 
are shown. The upper small figure shows the corresponding 
reduced of the ht. 


reflections of the FSI and the NE cusp in the pA invariant 
mass distributions contribute to the deviations. 

At the KE threshold no cusp-like structure is detected, 
as shown in fig. [T31 Here the KE threshold region is com¬ 
pared to the NE threshold region. The effect of the pE 
cusp is estimated to be (15 ± 5) % of the total cross sec¬ 
tion, which corresponds to (3 ± l)pb. From fig. [T3]it can 
be seen that the height of the NE cusp compared to the 
error bars of the measurement is about 10/1, therefore, 
the upper limit of a possible KE cusp - if it has a similar 
shape as the NE cusp - is in the order of 0.3 pb. 

The contribution of N* resonances decaying into KA 
will vary with the distance of the invariant mass KA to 
the resonance mass. As N* resonances affect the compo¬ 
sition of the angular momenta, this composition is tested 
for different ranges of the invariant mass KA. The angular 
distribution for the KA rest system given by the cosine 
of the Gottfried-Jackson angle is generated for six 

ranges of the KA invariant mass. The distributions and 
the corresponding ranges are shown in the appendix in 
fig. 1231 The angular distributions are fitted with a formula 
describing the differential cross section as a sum of Legen¬ 
dre polynomials by neglecting the interference terms. As 
the angular distributions do not have to be symmetric to 
cos'd = 0 - as it is the case of the center of mass system 
- the odd terms of the Legendre polynomials have to be 
included. 
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Fig. 14. The height of the cusp (green crosses) and the 
amount of the reflections at the cusp position (black crosses) 
for the different regions of cos ^. The blue dashed line shows 
the fit result of the height of the cusp of the total data set. As 
the errors of the fit parameter C/q are found to be asymmetric, 
the errors of the MIGRAD routine of the ROOT [37] program 
are shown. The upper small figure shows the corresponding 
reduced of fhe fit. 



invariant mass pA [GeV/c^] 


Fig. 15. The KE threshold region of the KA invariant mass 
is shown (top) in comparison to the NE threshold region of 
the pA invariant mass (bottom). The thresholds are indicated 
by the vertical lines. In contrast to the strong and sharp en¬ 
hancement at the NE threshold in the lower figure, no sharp 
structure is seen at the KE threshold in the upper figure. 


d " 

dcoJgXA = ^ (12) 

C is a normalization constant, for these fits the highest 
order n = 4 is chosen. The dependence of the Legendre 
coefficients Oi on the KA invariant mass slices is shown in 

fig. US) 

As the KA invariant mass is proportional to the K and 
A cm momentum the rising momenta of the KA invari¬ 
ant mass slices in addition influences the angular distri¬ 
butions, as it is shown in section o for the cm angular 
distributions. Therefore, an influence of N* resonances on 
the variation of the fit coefficients oq .. 04 with the in¬ 
variant mass can only be visible in the comparison with 
the corresponding fit coefficients calculated in the same 
manner for the invariant masses of pK and pA, where no 
resonances are expected. The GJ angular distributions for 
these invariant masses are shown in the appendix (figs. 1261 
E71) and the variation of the coefficients in figs. (TT] [TSl 



Fig. 16. Fit results of the coefficients of equation [12] obtained 
from the GJ angular distributions cos for six bins of the 

KA invariant mass with equal widths (lower figure). The coef¬ 
ficient uq is close to 1, it is given in fig. 1251 The upper figure 
shows the corresponding reduced of the fit. The dashed lines 
are generated with a second degree polynomial fit and have the 
only purpose to guide the eye. 

No specific behavior of the angular momentum com¬ 
position of the KA invariant mass is detected, which could 
be related to nucleon resonances. Therefore, this distribu¬ 
tion seems to be the wrong tool to detect the influence 
of nucleon resonances. The invariant mass spectra, which 
contain the A, which is the heaviest particle of the reac¬ 
tion in study, show both a nearly linear increase of the 02 
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Fig. 17. Fit results of the coefficients of equation ll2l obtained 
from the G J angular distributions cos 

Kp invariant mass with equal widths (lower figure). The coeffi¬ 
cient ao is close to 1, it is given in fig. 1261 The upper small figure 
shows the corresponding reduced of the fit. The dashed lines 
are generated with a second degree polynomial fit and have the 
only purpose to guide the eye. 


•o 

0) 2 
3 



coefficient with the invariant mass and a nearly quadratic 
increase of the 04 coefficient. For all three invariant masses 
the odd coefficients are close to zero, thus indicating that 
the angular distributions are almost symmetric. 


6 Summary 

More than 200 000 events of the reaction pp ^ pKA with a 
beam momentum of 2.95 GeV/c are analyzed. The angular 
distributions of the final state particles in the cm system 
exhibit different characteristics. They are analyzed in de¬ 
pendence on cm momentum by a fit of Legendre polynomi¬ 
als. While the Legendre coefficient 02 of the A distribution 
rises with the third power of the cm momentum, the co¬ 
efficient 02 of the proton distribution rises linear with the 
cm momentum and is approximately constant for values 
larger than 0.35GeV/c. The angular distribution of the 
kaon is nearly constant for all cm momenta. 

The threshold effect in the pA invariant mass spectrum 
at the NE threshold can be described by a Flatte distribu¬ 
tion by assuming a virtual state of pA with a mass above 
the threshold. It is shown by cuts on the pK helicity an¬ 
gle, that the height and the structure of the cusp is nearly 
constant along the KA invariant mass. Therefore, the in¬ 
tensity variations are most likely due to N* resonances 
decaying into KA. The height of the cusp is isotropic in 
dependence on the kaon cm angle and in dependence on 
the Gottfried-Jackson angle in the pA rest frame. 

In the KA invariant mass no deviations at the KS 
threshold are found. The upper-limit of a possible thresh¬ 
old effect is estimated to be 0.3 pb. The examination of 
the dependence of the Gottfried-Jackson angle in the KA 
frame on the KA invariant mass compared to the ones of 
the pA and Kp frames exhibits no peculiarities. This indi¬ 
cates that this examination may not be sensitive to effects 
of N* resonances. 
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Fig. 18. Fit results of the coefficients of equation ll2l obtained 
from the GJ angular distributions cos for six bins of the 
pA invariant mass with equal widths (lower figure). The coeffi¬ 
cient ao is close to 1, it is given in fig. 1271 The upper small figure 
shows the corresponding reduced of the fit. The dashed lines 
are generated with a second degree polynomial fit and have the 
only purpose to guide the eye. 
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Appendix 

In the appendix the individual fits of the observables, 
which are separated into bins of a second observable, are 
shown. Inside each figure the fit results and the range of 
the second observable are given. 
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Fig. 19. Proton angular distributions for different ranges of the proton center of mass momentum. The green curve shows 
the fit result. The contributions of the Legendre polynomials Pq, P 2 , and P 4 , weighted with the corresponding coefficient are 
plotted in black, red, and blue. 
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Fig. 20. Kaon angular distributions for different ranges of the kaon center of mass momentum. The green curve shows the fit 
result. The contributions of the Legendre polynomials Po, P 2 , and P4, weighted with the corresponding coefficient are plotted 
in black, red, and blue. The last figure is missing, as the maximum kaon momentum is less than 0.44GeV/c. 
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Fig. 21. A angular distributions for different ranges of the A center of mass momentum. The green curve shows the fit result. 
The contributions of the Legendre polynomials Pq, P2, and P4, weighted with the corresponding coefficient are plotted in black, 
red, and blue. 
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Fig. 22. The pA invariant mass distributions are shown for different slices of the Dalitz plot. The slices are generated with 
cuts in the helicity angle cosi?p^^ as given in the individual plots. The spectra are normalized to 1/6 of the total data set which 
is shown in fig. HOI The data are shown as points with error bars, the red line is the total fit. The green line is the contribution 
of the Flatte distribution, where only the height (CTo) is varied. The blue line is the FSI part with all parameters fixed. The 
black line is the part of the reflections with freely varying parameters of the second order polynomial. 
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Fig. 23. The pA invariant mass distributions are shown for different slices of the kaon cm scattering angle as given in the 
individual plots. The spectra are normalized to 1/6 of the total data set which is shown in fig. HOI The data are shown as points 
with error bars, the red line is the total fit. The green line is the contribution of the Flatte distribution, where only the height 
(C7~b) is varied. The blue fine is the FSI part with all parameters fixed. The black line is the part of the reflections with freely 
varying parameters of the second order polynomial. 















































20 


The COSY-TOF Collaboration: High statistics measurement of the reaction pp ^ pK^A 








Fig. 24. The pA invariant mass distributions are shown for different slices of the proton - beam angle in the GJ pA system 
as given in the individual plots. The spectra are normalized to 1/6 of the total data set which is shown in fig. HOI The data 
are shown as points with error bars, the red line is the total fit. The green line is the contribution of the Flatte distribution, 
where only the height (C/o) is varied. The blue line is the FSI part with all parameters fixed. The black line is the part of the 
reflections with freely varying parameters of the second order polynomial. 
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Fig. 25. The GJ angular distribution cos is shown for six slices of the KA invariant mass, the invariant mass ranges of 

the slices are given in the figures. The data, which are corrected to the acceptance and reconstruction efficiency, are shown by 
black error bars and the fit result according to equation [12] is shown by the red line. 
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Fig. 26. The GJ angular distribution cos is shown for six slices of the Kp invariant mass, the invariant mass ranges of 

the slices are given in the figures. The data, which are corrected to the acceptance and reconstruction efficiency, are shown by 
black error bars and the fit result according to equation [12] is shown by the red line. 
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Fig. 27. The GJ angular distribution cos shown for six slices of the pA invariant mass, the invariant mass ranges of 

the slices are given in the figures. The data, which are corrected to the acceptance and reconstruction efficiency, are shown by 
black error bars and the fit result according to equation [12] is shown by the red line. 
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